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Perturbative QCD-based models of parton production and equilibration in ultrarela- 
tivistic heavy ion collisions are reviewed with an emphasis on the treatment of quantum 
interference effects. Uncertainties in the initial parton production and their effects on later 
parton equilibration are considered. Probes of early parton dynamics are also discussed. 



1. INTRODUCTION 



In the last few years, there have been enormous new interests in a perturbative QCD- 
inspired description of ultrarelativistic heavy ion collisions based on the parton model [ |l| . 
^ ■ The main argument for such a treatment of high-energy heavy ion collisions lies in the the 
lyQ ', nuclear structure at different scales. When the transverse momentum transfer involved in 
^ I each nucleon-nucleon collision is small, pt ^ Aqcd, effective models based on,e.5^., meson- 
Q I exchange and resonance formation are sufficient to describe multiple interactions between 
' hadrons, in which parton structure of the hadrons cannot yet be resolved. Those coherent 



(with respect to partons inside a hadron) interactions lead to collective behaviors in low 
and intermediate energy heavy ion collisions as first observed in Bevalac experiments [ 0] 
and recently in experiments at the AGS energy [ |]. However, when becomes large 
enough to resolve individual partons inside a nucleon, the dynamics is best described on 
' the parton level via perturbative QCD (pQCD). Though hard parton interactions occur 
at CERN-SPS energies {^/s ^ 20 AGeV), they play a negligible role in the global features 
^ . of heavy ion collisions. However, at collider energies {^/s > 100 AGeV) the importance of 
hard or semihard parton scatterings is clearly seen in high-energy pp and pp collisions [ 
They are therefore also expected to be dominant in heavy ion collisions at RHIC and 
LHC energies [ |, . 

Hard or semihard interactions happen in a very short time scale and they generally 
break color coherence inside an individual nucleon. After the fast beam partons pass 
through each other and leave the central region, a dense partonic system will be left 
behind which is not immediately in thermal and chemical equilibrium. Partons inside 
such a system will then further interact with each other and equilibration will eventually 
be established if the interaction is frequent enough among the sufficiently large number 
of initially produced partons. Due to the asymptotic behavior of QCD, production rates 
of hard and semihard partons are calculable via pQCD during the initial stage of heavy 
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ion collisions. The color screening mechanism in the initially produced dense partonic 
system may make it also possible to use pQCD to investigate the thermal and chemical 
equilibration of the system. In this talk, I will give a brief review of parton production in 
high energy heavy ion collisions. Then I will concentrate on recent developments in early 
parton dynamics. I will also discuss the uncertainties in the initial condition of parton 
equilibration and their consequences in the formation of a quark-gluon plasma. 

2. INITIAL PARTON PRODUCTION 

In a parton model, a nucleus can be considered as an assemble of interacting partons. 
Because of vacuum fluctuations, a valence quark can radiate many off-shell gluons, quarks 
and antiquarks. These partons, if undisturbed, will reassemble back to the valence quark. 
However, a hard scattering can break the coherence and bring these virtual partons on 
shell. The parton number density fa/A{x, Q"^) thus depends on the resolution of the probe 
or the momentum transfer of the collision. Similarly, a hard scattering can also lead to 
final state radiations off the scattered partons which can be described by fragmentation 
functions Da/biZyQ"^)- In the leading logarithmic approximation in an axial gauge, the 
interference between multiple radiations and between initial and final state radiations can 
be neglected. The amplitude for successive radiations then has a simple ladder structure 
and both the parton distribution functions fa/Ai^iQ'^) ci^nd the fragmentation functions 
follow the Altarelli-Parisi evolution equations [ 0] with respect to Q^. 

Using the parton distribution functions fa/A{x, Q^), parton scattering rates in nucleus- 
nucleus collisions can be written as 



where daab is the cross section for parton-parton scatterings, yi and y2 are the rapidities 
of the scattered partons, xi and X2 are the light-cone momentum fractions carried by 
the initial partons, and the summation runs over all parton species. The factor K 2 
accounts for next-to-leading order effects. The parton distribution density of a nucleus is, 

fa/A{.X, r) = tA{j)Sa/A{.X, v)fa/N{.X, Q^), (2) 

where ^^(r) is the thickness of the nucleus which is normalized to / (i^rtA(r) = ^, 
fa/N{x,Q'^) is the parton structure function of a nucleon and Sa/A{x,r) accounts for the 
nuclear modification of parton distributions. 

Many models have been developed to simulate parton production based on Eq. ([l|) ^ 
|T0|, ^ . In particular, HIJING and PCM have applied the techniques of renormalization 
group improved perturbative QCD, initially developed to study jet formation in e^e~ and 
pp collisions [ [1^, |I3] , to simulate initial and final state radiations associated with hard 
or semihard parton scatterings. Shown in Fig. || are the rapidity densities of produced 
partons given by HIJING and PCM [H. Both models predict a large number of produced 



partons which will form a dense partonic gas as the initial condition immediately after 
the overlap of the two colliding nuclei with formation time, l/pr ~ 0.2 fm/c. To ensure 
the credibility of these predictions, all models have been checked against the existing data 
of pp and pp collisions. It has been demonstrated that many aspects of multi-particle 
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Figure 1. Rapidity distributions of produced partons given by HIJING and PCM calcu- 
lations for central Au + Au collisions at RHIC energy 



production in pp and pp collisions can be accounted for by the onset of semihard parton 
scatterings in high energy hadronic collisions. 

Even though existing pp and pp data can provide many constraints on these pQCD mo- 
tivated phenomenological models, they still cannot uniquely determine all the parameters, 
giving rise to many uncertainties in the initial parton production: 

(l)To regulate the infrared divergences in the QCD cross section of minijet production, 
a pt cut-off po has to be introduced. This cut-off and the corresponding soft inclusive 
cross section cTsoft for interactions with smaller pt < Po are constrained by the energy 
dependence of the total pp and pp cross sections. These two parameters cannot be de- 



termined uniquely from the phenomenology of currently measured observables [15]. To 
reduce the uncertainties, measurement of two-particle correlation functions in azimuthal 
angle in the transverse plane was proposed [ |16|] to further constrain po and cxsoft- 

(2) Since most of produced partons are gluons, their number is sensitive to the gluon 
distribution function fg{x) at small x. So far there is still no precise measurement of 
the gluon distribution function in the small x region where most of the minijets originate 
at RHIC and LHC energies. Though there have been theoretical attempts to study 



gluon distribution in a nucleus [ there is no direct experimental measurement of 

nuclear shadowing of the gluon distribution. Future experiments at HERA might provide 
information on the gluon distribution. In addition systematical measurements of pp, pA 
and AA interactions at RHIC can be used to study gluon shadowing [ 1^ . 

(3)The treatment of primary-secondary parton scatterings can also lead to differences in 
the initial parton production. Because of formation time effect which I will discuss later. 
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Values of the relevant parameters characterizing the parton gas at time Ti, 
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most of the produced partons have no time to rescatter again with another oncoming 
beam parton before they pass through. Proper inclusion of such formation time effect 



will suppress primary-secondary parton scatterings [20, 21 1. Systematic studies of pA 
collisions at RHIC will provide a quantitative estimate of this formation time effect. 

(4) Finally, parton production from soft processes is also very model dependent. In 
parton cascade model, soft interactions are treated as elastic parton scatterings with a 
regularized distribution at small pr- Such soft interactions with large cross section cXgQft 
can lead to many produced partons, especially when the long interaction time is not taken 
into account. In models with a string phenomenology, it is argued that soft interactions, 
though not calculable in pQCD, produce coherent gluons which form strings between 
leading partons. These strings can materialize into soft partons but at rather longer time 
scale of more than 1 fm/c. Color screening by the hard partons produced earlier can also 
reduce the soft parton production from the coherent color field [ |2^ . 

3. EARLY PARTON DYNAMICS 

The formation time of parton production can be estimated via uncertainty principle. 
Thus model calculations can also give the space-time history of parton production. Not 



surprisingly, it is found [ ^ that most of partons are produced within 0.4 fm/c after the 
complete overlap of two colliding nuclei. The parton system will then undergo further 
interactions and free-streaming. Neglecting parton rescatterings in this period of time, 
the kinematic separation of partons with different rapidities in a cell establishes local 
momentum isotropy at the time of the order of Tiso [ 0, Given this initial time Tiso, 
the parton rapidity density dN/dy and the average transverse momentum (/ct), the initial 
parton number and energy density can be estimated via Bjorken formula. 



dN 4,, , . , 

"° = ^^^' ^o = -o-(M- (3) 



Assuming further that the phase distributions have a factorized form, fi{k) = Aj/eq(fc,T), 
one can then estimate the effective initial temperature and fugacities. Listed in Table 1 
are the initial conditions for central Au + Au collisions at RHIC and LHC collider energies 
from HIJING calculations. 



One can observe that the initial parton gas is rather hot reflecting the large average 
transverse momentum. However, the parton gas is still undersaturated as compared to 
the ideal gas at the same temperature. The gas is also dominated by gluons with its quark 
content far below the chemical equilibrium value, confirming the hot glue scenario [ 

How this parton gas evolves toward equilibrium is under intense investigation in the 
last few years. Tremendous progresses have been made with different approaches. Parton 
cascade model looks particularly interesting since it can provide a space-time picture of 
the evolution toward a thermalized quark-gluon plasma. Because of the classical nature 
of the cascade model, however, it is rather difficult to simulate many subtle quantum 
interference effects which are expected to be very important. For example, M. Gyulassy 
and I have studied the criteria for a semiclassical treatment of multiple parton scatterings. 
We found [ ^ that the billiard ball picture of multiple scatterings in a cascade model is 
valid only when the momentum transfer of each scattering is large enough to resolve the 
spatial separation of nearby scatterers which may not be true in the initial scatterings of 
the beam partons. Another important problem is the color interference in QCD and how 
it can be incorporated into semiclassical simulations. 

The importance of color interference can be best demonstrated by comparing induced 
radiation in QED and QCD by a single scattering. In QED case, the radiation amplitude 
can be separated into initial state and final state radiation in the soft radiation limit, 

,J. = »e(p-fW), (4) 

where k and e are the momentum and polarization of the radiated photon, respectively. 
Pi and pf are the initial and final momentum of the charged particle. In the high energy 
limit, the two terms in TZi cancel to the order 0{l/^/s) in the central region due to 
the destructive interference between initial and final state radiations. This gives rise to 
a valley in the photon's rapidity distribution as illustrated in Fig. |^(a). However, the 
cancellation will not be complete in QCD because gluons from initial and final state 
radiations have different colors due to the color exchange in the scattering. Furthermore, 
a gluon can also be emitted from the gluon propagator via the three gluon vertex. The 



resultant gluon spectrum [ |26 

dyd?k^ ~ 7r2 fci(qx-kx)2' ^ ' 

has a plateau in the central rapidity region, as shown in Fig. 0(b). This is very different 
from the QED case. 

Similarly, destructive interference between induced radiations by neighboring scatter- 
ings is also important. The radiation amplitude induced by multiple scatterings in QED, 
for example, is, 

^™-e|e-(i^-i^). (a) 

In the eikonal limit ( particle propagating along a straight line), the phase difference 
between any two terms can be written as, 

k.{x,- X,) = U,/r{k)- r{k) = ^j^^^^ - ^, (7) 




Figure 2. Illustrations of radiation spectrum in (a) QED and (b) QCD induced by a single 
scattering. 



where Lij is the distance (time) between scattering i and j, and r(/c) is the effective 
formation time. In the Bethe-Heitler limit Lij 3> t(/c), the total intensity of induced 
radiation is simply additive in the number of scatterings. However, when L^j <^ T{k), the 
final state radiation of a scattering completely cancels with the initial state radiation from 
the previous one. This so-called Landau-Pomeranchuk-Migdal (LPM) effect can also be 
understood in terms of uncertainty principle: A radiated photon ( or gluon in QCD) must 
travel at least one wavelength l/k± in the transverse direction in order to be separated 
from its parent particle. The radiation will be suppressed if the parent particle suffers 
another scattering within time (distance) r(/c), leading to the LPM effect. 

Analysis in QCD gives similar results, although one has to take into account different 
color structures associated with different scatterings. In Ref. [ |2^, M. Gyulassy and 
I considered induced radiation by multiple scatterings in the soft radiation limit and 
neglected the situation where a radiated gluon can suffer further interactions. We obtained 
a suppressed gluon spectrum with effective formation time, 

^QCd(A:) = ^r{k). (8) 



Baier et al. [ recently argued that those diagrams which were neglected in Ref. [ ^ 
might be important especially when a radiated real gluon travels with the beam particle for 
a long time and then is knocked off later. This argument will have important consequences 
in the energy loss of a fast parton propagating in a QCD medium [ However, I will 
use the above analysis in terms of effective formation time to discuss induced radiation 
in the context of parton equilibration. 



4. PARTON EQUILIBRATION RATE 



I would like to emphasize that the above interference effects happen on the matrix 
elements level. In particular, the LPM effect involves the interference between initial 



and final state radiations from different scatterings. A detailed analysis [ 25, ^ can 
show that the interference actually happens between two amplitudes in which the beam 
parton has completely different virtualities, i.e., time-like in the final state radiation of one 
scattering and space-like in the initial state radiation of the previous one. This imposes 
a great difficulty for a proper treatment of the quantum interference in a parton cascade 
model in which a parton must always remain time-like between two scatterings [ A 
possible remedy for this is to consider both initial and final state radiation associated with 
each scattering together and include the interference effect by using a modified radiation 
spectrum. The effective spectrum should suppress soft gluon radiation whose formation 
time is larger than the mean-free-path of parton scatterings. To demonstrate this, let us 



consider parton equilibration in the form of rate equations [ . The analysis of QCD 
LPM effect [ ^ ^ has been done for a fast parton traveling inside a parton plasma. I 
will apply the results to radiations off thermal partons whose average energy is about T, 
since we expect the same physics to happen. 

We consider only the dominant process gg ggg in the leading order. In order to per- 
mit the approach to chemical equilibrium, the reverse process, i.e., gluon absorption, has 
to be included as well, which is easily achieved making use of detailed balance. Radiative 
processes involving quarks have substantially smaller cross sections in pQCD, and quarks 
are considerably less abundant than gluons in the initial phase of the chemical evolution 
of the parton gas. Here we are interested in understanding the basic mechanism under- 
lying the formation of a chemically equilibrated quark-gluon plasma, and the essential 
time-scale. We hence restrict our considerations to the dominant reaction mechanism for 
the equilibration of each parton flavor. These are just four processes: 

99 ^ 999, 99 ^ (9) 

Other scattering processes ensure the maintenance of thermal equilibrium [gg ^ gg, gq ^ 
gq, etc.) or yield corrections to the dominant reaction rates {gq ^ qgg, etc.). 

Restricting to the above reactions and assuming that elastic parton scatterings are suf- 
ficiently rapid to maintain local thermal equilibrium, the evolution of the parton densities 
is governed by the master equations [ ^ : 

a,Xn,u^) = (11) 

where hi = Ui/Xi denote the densities with unit fugacities. A, = 1, (T3 and (T2 are thermally 
averaged, velocity weighted cross sections, 

^3 = W{99 999)v), 02 = Wigg qq)v). (12) 



We have also assumed detailed balance and a baryon symmetric matter, = Uq. If we 



neglect effects of viscosity due to elastic scattering [ , we then have the hydrodynamic 
equation 

^ + ^ = 0. (13) 

We further assume the ultrarelativistic equation of motion, e = 3P. We can then solve 
the hydrodynamic equation, 

[Ag + — A,]=^/^rV = const. , (14) 

and rewrite the rate equation in terms of time evolution of the parameters T(r), Ag(r) 
and Ag(r), 

^ + 3^ + ^ = i?3(l-A,)-2i?2(l-|) (15) 
where the density weighted reaction rates R3 and R2 are defined as 

^3 = ^CTsng, R2 = \(^2ng. (17) 

Notice that for a fully equilibrated system (A^ = A, = 1), Eq. ( p!4D corresponds to the 
Bjorken solution, T{t)/Tq = (tq/t)^/^. 

To take into account of the LPM effect in the calculation of the reaction rate for 
99 ~^ 999i W6 simply impose the LPM suppression of the gluon radiation whose effective 
formation time tqcd is much longer than the mean-free-path A/ of multiple scatterings 
to each gg ggg process. In the mean time, the LPM effect also regularizes the infrared 
divergency associated with QCD radiation. However, still contains infrared singularities 
in the gluon propagators. For an equilibrium system one can in principle apply the 



resummation technique developed by Braaten and Pisarski [ |SD| to regularize the electric 
part of the propagators, though the magnetic sector still has to be determined by an 
unknown magnetic screening mass. Since we are dealing with a nonequilibrium system, 
Braaten and Pisarski's resummation may not be well defined. As an approximation, we 
will use the Debye screening mass [ ^ , 



1^1 = -^ r kf{k)dk = AirasT^Xg, (18) 
vr^ JO 

to regularize all singularities in the radiative cross section. 

To further simplify the calculation we approximate the LPM suppression factor in Ref. [ 
^] by a 6'-function, 6'(A/ — tqcd)- The modified differential cross section for gg —>■ ggg 
is then, 



(19) 



Figure 3. The scaled 
gluon production 
rate R^/T (solid line) 
for gg ggg and 
the quark production 
rate R2IT (dashed 
line) for gg qq 
as functions of the 
gluon fugacity Xg for 
a, = 0.3. 
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where the second step-function accounts for energy conservation, and s = 18T^ is the 
average squared center-of-mass energy of two gluons in a thermal gas. The regularized 
gluon density distribution induced by a single scattering is, 

dug _ Ca^s ql ^20) 



Similarly, the regularized small angle gg scattering cross is, 
da'J 9 27ral 

dql ^{ql + ^^lr ^ ' 

The mean-free-path for elastic scatterings is then, 

A7^ = riq / dq\—^ = -aiUsT -— , (22) 

which depends very weekly on the gluon fugacity A^. 

We can evaluate the integration numerically and find out the dependence of R3/T on 
the gluon fugacity A^ as shown in Fig. 0, for a coupling constant = 0.3. The gluon 
production rate increases with Xg and then saturates when the system is in equilibrium. 
Overall, the gluon equilibration rate is suppressed due to the inclusion of the LPM effect 
on induced radiation. 

The calculation of the quark equilibration rate R2 for gg — > qq is more straightforward. 
Estimate in Ref. [ ^ gives, 

R2 = ^cr^Ug ^ OMNf a%T ln{5. 5 /Xg). (23) 

The dashed line in Fig. ^ shows the normalized rate R2/T for Nf = 2.5, taking into account 
the reduced phase space of strange quarks at moderate temperatures, as a function of the 
gluon fugacity. 
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Figure 4. Time evolution of the temperature T and the fugacities A^, and \g of gluons and 
quarks in the parton plasma for central Au + Au collisions at the (a) RHIC and (b) LHC 
energies. The initial values for T, and Xq are determined from HIJING simulations 
and are listed in Table ^ 



5. EVOLUTION OF THE PARTON PLASMA 

With the parton equilibration rates which in turn depend on the gluon fugacity, we 
can solve the master equations self-consistently and obtain the time evolution of the 
temperature and the fugacities. Shown in Fig. ^ is the time dependence of T, \g, and A^ 
for initial conditions listed in Table |I] at the RHIC and LHC energies. We find that the 
parton gas cools considerably faster than predicted by Bjorken's scaling solution (T^r = 
const.) shown as dotted lines, because the production of additional partons approaching 
the chemical equilibrium state consumes an appreciable amount of energy. The accelerated 
cooling, in turn, slows down the chemical equilibration process, which is more apparent 
at the RHIC than at LHC energies. Therefore, the parton system can hardly reach its 
equilibrium state before the effective temperature drops below Tc ~ 200 MeV in a short 
period of time of 1-2 fm/c at the RHIC energy. At the LHC energy, however, the parton 
gas becomes very close to its equilibrium and the plasma may exist in a deconfined phase 
for as long as 4-5 fm/c. Another important observation is that quarks never reach chemical 
equilibrium at both energies. This is partially due to the small initial quark fugacity and 
partially due to the small quark equilibration rate. 

In another recent development, E. Shuryak and L. Xiong [ ^ used the complete ma- 
trix elements of gg gg + {n — l)g which include complete interference to calculate the 
parton equilibration rate. This is essentially a step beyond leading logarithmic approxi- 
mation used in most of the calculations. They found that inclusion of these higher order 
corrections not only increases the initial parton production but also the parton equilibra- 




Figure 5. The same as in Fig. Ka), except that the initial parton densities are 4 times 
higher than given in Table with the same (ordinary lines), or reduced initial temperature, 
To = 0.4 GeV (lines with stars) 



tion rate. However, it is expected that if the formation time and virtual corrections (or 
unitarity) are included in this multiple gluon radiation processes, the effect could become 
smaller. 

As I have pointed out, there are still a lot of uncertainties in the initial parton produc- 
tion. We can estimate the effect of the uncertainties in the initial conditions on the parton 
gas evolution by multiplying the initial energy and parton number densities at the RHIC 
energy by a factor of 4. This will result in the initial fugacities, A° = 0.2 and Ajj = 0.024. 
With these high initial densities, the parton plasma can evolve into a nearly equilibrated 
gluon gas as shown in Fig. ^. The deconfined phase will also last longer for about 4 fm/c. 
The system is still dominated by gluons with few quarks and antiquarks as compared to 
a chemically equilibrated system. If the uncertainties in the initial conditions are caused 
by soft parton production from the color mean field, the initial effective temperature will 
decrease. Therefore, we can alternatively increase the initial parton density by a factor 
of 4 and decrease Tq to 0.4 GeV at the same time. This leads to higher initial fugacities, 
\g = 0.52 and = 0.083. As shown in Fig. | by the curves with stars, this system 
evolves faster toward equilibrium, however, with shorter life-time in the deconfined phase 
due to the reduced initial temperature. 

We thus can conclude that perturbative parton production and scatterings are very 
likely to produce a quark-gluon plasma ( or more specifically a gluon plasma) in ultra- 
relativistic heavy ion collisions at the RHIC and LHC energies. However, uncertainties in 
the initial conditions have to be carefully examed in order to give a definite prediction of 
the quark-gluon plasma at the RHIC energy. 



6. PROBES OF EARLY PARTON DYNAMICS 



While some of the uncertainties may be reduced by theoretical studies of the parton 
distribution in nuclei and initial parton production, most of them have to be resolved 
experimentally. There are many processes the measurements of which can provide exper- 
imental probes of the early parton dynamics. For example, suppression of large pt jet or 



jet quenching [ |T^, |33|, 0, |3^ due to the final state interaction of large pt parton with 
the parton gas can be used to study parton equilibration in the early stage. Furthermore, 
the energy loss of the jets and the alcoplanarity [ ^ |3^ can also be used to study the 
thermal dynamical properties of the parton gas . 

Among these hard probes, electromagnetical signals like thermal photons and dileptons 
are considered more direct since they can escape the dense matter without further inter- 
actions. They can thus reveal the dynamics of initial parton production and equilibration. 
In addition, open charm production can also be considered as a direct probe since charm 
quarks cannot be easily produced during the mixed and hadronic phases of the dense 
matter due to their heavy masses. To leading order in pQCD, dilepton production is 
dominated by qq direct photon by q{q)g q{q)l and open charm by gg cc 

processes. As pointed out by Strickland [ ^ , measurements of these direct thermal sig- 
nals can tell us the relative ratio of quark and gluon number densities in the early stage 
of parton equilibration. 

To demonstrate the sensitivity of these direct probes to the initial condition of the 
parton evolution, let us consider open charm production as an example. Charm produc- 
tion can be divided into three different contributions in the history of the evolution of 
the parton system: (1) initial production during the overlapping period ; (2) pre-thermal 
production from secondary parton scatterings during the thermalization, r < Tiso; (3) 
and thermal production during the parton equilibration, r > riso, in the expanding sys- 
tem. The initial charm production can be calculated similarly to minijet production [cf. 
Eq. ([l|)]. For thermal and pre-thermal production, the differential rate is[ 

73 



a-^p 16(27r)** J u^ uj2 E2 

(24) 



l9hfAh)fgik2)\M,,^,-f + g'jqik,)f,ik2)\M,,^M' 

given the phase-space density of the equihbrating partons, fi{k), where gG=16, gq = 6A/, 
are the degeneracy factors for gluons and quarks (antiquarks) respectively, |A^gg^ccP, 
\A4gg^cc\'^ are the averaged matrix elements for gg cc and qq cc processes. Due to 
small charm density, the Pauli blocking of the final charm quarks can be neglected. The 
corresponding charm spectrum is 

dNc d2 /■ - - ^d^^ 



t:R\ / drjdTE——. (25) 



dyd'^p± J d^p 

For pre-thermal charm production, one can use the momentum spectrum of initially 
produced partons and model the phase distribution by introducing a correlation between 
momentum and space-time. This correlation was shown to be important in the calculation 
of pre-thermal charm production [ ^]. For thermal production, one can use the time 




Figure 6. The p± distribution of the initial (sohd), pre-thermal (dot-dashed) and ther- 
mal (dashed) charm production for initial parton densities 4 times higher than HIJING 
estimate given in |l| but with the same (ordinary lines), or reduced initial temperature, 
To = 0.4 GeV (lines with stars). 



dependence of temperature and fugacities as given by the parton evolution equations. 
Shown in Fig. ^ are the charm production rates as functions of pt with the same initial 
conditions in Fig. |^. Note that with the increased initial parton density, open charm 
production from pre-thermal and thermal stage is as important as the initial production. 
Thermal production is more sensitive to the variation of the initial temperature than the 
initial fugacities. Therefore, by measuring charm enhancement, we can probe the initial 
parton phase-space distribution, initial temperature and equilibration time. 

7. CONCLUSIONS 

In this talk I reviewed the pQCD-based models of ultrarelativistic heavy ion collisions. 
In this framework the reaction dynamics can be described by perturbative parton scat- 
terings and radiations. Using model estimate of initial parton production, it is found 
that there are enormous number of partons produced during the overlap period of two 
colliding nuclei. Thereafter, local isotropy in momentum distribution might be reached. 
Parton proliferation through induced radiation and parton fusion further drive the parton 
system toward a fully equilibrated parton plasma. The quarks always lag behind gluons 
in reaching their equilibrium. Due to the energy consumption by parton production, the 
parton gas cools down faster leading to a reduced plasma life time of 4 - 6 fm/c. 

Throughout the production and evolution of the partonic system, interference effects 
play an important role in multiple collisions, the correct implementation of which will be 
a great challenge to any classical cascade models. Especially, the destructive interference 
among different amplitudes of gluon radiation induced by multiple scatterings suppresses 



soft gluons whose formation time is larger than the mean-free-path of parton scatterings 
inside a QCD medium. A detailed analysis of the radiation amplitudes reveals the under- 
lying physics which contradicts the intuitive picture of a classical cascade model. One way 
to incorporate the LPM effect in the parton interaction simulations is to consider both 
the initial and final state radiation together for each scattering and impose the formation 
time requirement, tqcd < ^f^ for the integration over the phase-space of the radiated 
gluons. This will lead to a reduced parton equilibration rate. 

I also discussed the uncertainties in the initial conditions and the experimental probes 
of the early parton evolution. All in all, these uncertainties arise from our ignorance 
of the nonperturbative physics and our inability to calculate the soft processes in the 
framework of QCD. In the pQCD-based models I have reviewed here, the uncertainties 
really lie in the cut-off, po, which supposes to separate nonperturbative soft interactions 
from perturbative hard processes. Since soft and hard physics do not have a definite 
boundary, the resultant parton production from hard or semihard is very sensitive to the 
cut-off. The accompanying soft parton production is not known in this model and may 
only be estimated by simple models like the color flux-tube model [ ^ . There have been 
recent developments in the field theory of particle production from mean-fields [ ^]. 
The chaotic behavior of nonabelian gauge fields may be intimately related to multiple 
parton production and fast gluon equilibration [ Q . Such a field theoretical approach to 
particle production could be the ultimate and consistent way to address the production 
and formation of a parton plasma in heavy ion collisions. 

In conclusion, parton equilibration is an exciting new subject with many unsolved 
problems. I hope eventually, not far from now, we can explicitly explain whether and how 
a quark-gluon plasma is formed in a heavy ion collisions. 
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